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The weak equivalence principle of gravity is examined at the quantum level in two ways. First, 
the position detection probabilities of particles described by a non-Gaussian wave-packet projected 
upwards against gravity around the classical turning point and also around the point of initial 
projection are calculated. These probabilities exhibit mass-dependence at both these points, thereby 
c 2 ^ , reflecting the quantum violation of the weak equivalence principle. Secondly, the mean arrival time 
04 ' of freely falling particles is calculated using the quantum probability current, which also turns 
fj ' out to be mass dependent. Such a mass-dependence is shown to be enhanced by increasing the 
^ , non-Gaussianity parameter of the wave packet, thus signifying a stronger violation of the weak 
f^ • equivalence principle through a greater departure from Gaussianity of the initial wave packet. The 
^^ ' mass-dependence of both the position detection probabilities and the mean arrival time vanish in 
C^^ , the limit of large mass. Thus, compatibility between the weak equivalence principle and quantum 
Cn . mechanics is recovered in the macroscopic limit of the latter. A selection of Bohm trajectories is 
exhibited to illustrate these features in the free fall case. 
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, ^ ' I. INTRODUCTION 

(N 
> 

lO , In the famous gedanken experiment conceived by Galileo, the universality of the ratio between the gravitational 

^~~l • and inertial masses has been studied with test bodies in free fall from the leaning tower of Pisa [l| . Since then several 
^*l [ tests of weak equivalence principle have been performed with classical test bodies such as very sensitive pendula or 
t~^ ' torsion balances establishing extremely accurate conformation of the equality of gravitational and inertial masses. 
^^ I The most familiar tests of the weak equivalence principle are experiments of the Eotvos-type [3|, which measure 
the gravitational acceleration of macroscopic objects. Thus Einstein's weak equivalence principle is established in 
the Newtonian picture of the gravitational force which is exactly proportional to the inertial mass. The validi_ty_of 

the principle even for quantum mechanical particles is also proved using gravity-induced interference experiments 
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The traditional equivalence principle is fundamentally both classical and local, and it is interesting to enquire 
j^ • how it is to be understood in quantum mechanics. Principle of equivalence states that inertial mass is equal to 
gravitational mass; mi = rrig = m. Another statement is that when all sufficiently small test bodies fall freely, they 
acquire an equal acceleration independent of their mass or constituent in a gravitational field. Quantum mechanically, 
it [5| can be said that " The results of experiments in an external potential comprising just a sufficiently weak, 
homogeneous gravitational field, as determined by the wave function, are independent of the mass of the system" . 
The last statement is also known as weak equivalence principle of quantum mechanics (WEQ). 
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The possibility of violation of weak equivalence principle in quantum mechanics is discussed in a number of papers, 
for instance using neutrino mass oscillations in a gravitational potential Q. The evidence of violation of WEQ 
can be shown both experimentally and theoretically. Experimental evidence exists in the interference phenomenon 
associated with the gravitational potential in neutron and atomic interferometry experiments (3|, |j| . Theoretically, 
for a particle bound in an external gravitational potential, it is seen that the radii, frequencies and binding energy 
depend on the mass of the bound particle [3-[al ■ 

The free fall of test particles in a uniform gravitational field in the case of quantum states with and without a 
classical analogue was re-examined by Viola and Onofrio[10j. Another quantum mechanical approach of the violation 
of WEQ was given by Davies [111 ] for a quantum particle in a uniform gravitational field using a model quantum clock 
[12| . A quantum particle moving in a gravitational field can penetrate the classically forbidden region beyond the 
classical turning point of the gravitational potential and the tunnelling depth depends on the mass of the particle. 
So, there is a small mass-dependent quantum delay in the return time representing a violation of weak equivalence 
principle. This violation is found within a distance of roughly one de Broglie wavelength from the turning point of 
the classical trajectory. 

In the recent past, Ali et al. [l3| have shown an explicit mass dependence of the position probabilities for quantum 
particles projected upwards against gravity around both the classical turning point and the point of initial projection 
using Gaussian wave packet. They also have shown an explicit mass dependence of the mean arrival time at an 
arbitrary detector location for a Gaussian wave packet under free fall. Thus both the position probabilities and the 
mean arrival time show the violation of WEQ for quantum particles represented by Gaussian wave packets. In the 
present paper, using a class of non-Gaussian wave packets involving a tunable parameter, we shall compute position 
detection probabilities around the classical turning point and the point of initial projection and also compute the 
mean arrival time for atomic and molecular mass particles represented by non-Gaussian wave packets to observe 
stronger violation of WEQ for quantum particles. To compute the mean arrival time, we consider the case when the 
particles are dropped from a height with zero initial velocity. 

Violation of WEQ arises as a consequence of the spread of wave packets the magnitude of which itself depends on 
the mass. In order to illustrate this effect we present an analysis in terms of Bohmian trajectories. In the Bohmian 



model [IJ] each individual particle is assumed to have a definite position, irrespective of any measurement. The 
pre-measured value of position is revealed when an actual measurement is done. Over an ensemble of particles having 
the same wave function ^, these ontological positions are distributed according to the probability density p = |^p 
where the wave function -0 evolves with time according to the Schrodinger equation and the equation of motion of any 
individual particle is determined by the guidance equation v = J / p^ where v is the Bohmian velocity of the particle 
and J is the probability current density. Solving the guidance equation one gets the trajectory of the particle. In the 
present work we describe the evolution of Bohmian trajectories in terms of the mass and the tunable parameter of 
the non-Gaussian wave packet. 

II. POSITION DETECTION PROBABILITY 



Let us consider an ensemble of quantum particles projected upwards against gravity with a given initial mean 
position and mean velocity. Our motivation is to look for greater violation of the equivalence principle compared to 
that achieved by a Gaussian wave packet. Hence, we formulate a one-dimensional non-Gaussian wave function which 
represents the initial state of each particle, given by 
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Here, a is the tuneable parameter varying from to 1. The normalisation factor N is given by 
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The salient features corresponding to the above wave packet are its asymmetry, its infinite tail, and its reduciblity to 
a Gaussian wave packet upon a continuous decrement of the parameter a to zero. The property of an infinite tail of 



the above wave function is not associated with non-Gaussian forms that are generated by truncating the Gaussian 
distribution. The asymmetry of the wave packet (due to the sine function) entails a difference between its mean and 
peak, as will be further evident in our following analysis. 

The initial group velocity is given by u = — -, where m is the mass of the particle and fco is the wave number. For 
a = 0, we get the Gaussian wave function. If ipi and ip2 represent the initial wave functions for particles 1 and 2 in 
the Schrodinger picture, then the mean initial conditions are 
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where (z)^ and (pz)^ represent the expectation values for position and momentum operators respectively. Here, we 
have restricted ourselves to a one-dimensional representation along the vertical z direction for simplicity. With the 
propagator 
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of a particle in the linear gravitational potential V = mgz and with the relation 
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one can get the Schrodinger time evolved wave function ip(z,t) at any subsequent time t as 
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where, St = fio ( 1 + » , ^ ^ ) ^^^ /^ ^ T^- The probability density is given by 
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FIG. 1: The variation of the width of the wave packet with a for different constant values of m (in a.m.u.). We take cro = 10 ^ 
cm, f = 1 sec and u = 10"^ cm s^^. 
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is the spreading of the Gaussian wave packet. The probability density is explicitly mass dependent. Here, the 
spreading of the wave packet is given by 



where. 
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which is also mass dependent. 

Fig. [1] gives the plot of a^G vs a. 

From this figure, we can see that for large masses, (T]\[q becomes almost constant with a. For a particular value of 
a, aNG decreases with increasing value of mass m. 



The expression for expectation value of z is given by 
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2 + a^(l~e~ 
Ehrenfest's theorem states that expectation value of quantum mechanical operators obey classical laws of motion. 
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FIG. 2: The variation of Zpeak and (z) with time t. We take (Tq = 10 
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50 a.m.u. and a = 0.5. 
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the instances of which are given below 



(12) 

It is possible to check that Ehrenfest's theorem is satisfied for our non-Gaussian wave packet. In the Gaussian limit 
(i.e., a — > 0), (z) = Zp^ak and the dynamics of (z) is like a classical point particle, as expected. 

In Eq. (jlip the first term is a-dependent. There is no m in this expression. So, (z) is a-dependcnt but mass 

independent. With the solutions obtained from the equation -'j-^ = 0, we can get the value of Zpe^fc for which |^p 
is maximum. But the equation can't be solved analytically. So, we have computed Zp^ak numerically. Zpeafc depends 
both on a and mass. In Table 1, it is shown numerically how Zpeak and (z) vary with a, for < a < 1. In Table 2, it 
is shown numerically how those two vary with mass. It is clear that Zpeak 7^ (z) and Zpeak increases with a and mass 
m but (z) increases with a, and remains constant for all masses. 

Fig. [2] shows the variation of Zpeak and (z) with time simultaneously. Dotted curve shows the motion of Zpeak 
whereas the continuous curve shows the motion of (z). As we take a^^ = 10"'' cm, the a-dependence of (z) and Zpg^ak 
becomes vanishingly small. So, this a-dependcncc can not be shown graphically. As a result, the curves for Zpeak and 
(z) coincide totally in this graph. 

Now, we want to observe the expression for (z) in the classical limit. If a full classical description is to emerge 
from quantum mechanics, we must be able to describe quantum systems in phase-space. The Wigner distribution 
function |15l . Ilq which is one of the phase-space distributions, can provide a re-expression of quantum mechanics in 
terms of classical concepts. 



The Wigner distribution function which is calculated from the time evolved wave function ^(z,t), is given by 
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and the marginals of which yield the correct quantum probabilities for position and momentum separately. 



By substituting the value of i^*{z -f- y,t) and ^-){z — y,t) using Eq.(5) we obtain the Wigner distribution function 
D^{z,p,t) for non-Gaussian wave packet. It satisfies the classical Liouville's equation [17| given by 
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TABLE I: The variation of (2) and Zp^ak in cm with a for mass tti = 10 a.m.u. t — 2 sec and cro = 0.1 cm. 



a Zpeak 

in cm 


in cm 


38.59999999999815 


38.599999999999910 


0.1 38.61540633397262 


38.611498366598200 


0.2 38.62925626009014 


38.622755653868810 


0.3 38.64090000000000 


38.633500000000000 


0.4 38.65026513575232 


38.643679695976970 


0.5 38.65786585282045 


38.652999876917650 


0.6 38.66400000000000 


38.661400000000000 


0.7 38.66910000000000 


38.668800000000000 


0.8 38.67337906376510 


38.675246198243485 


0.9 38.67696171641449 


38.680689435211890 


1.0 38.68002216630916 


38.685197660661515 



TABLE IL The variation of (z) and Zpeak in cm with mass m in a.m.u. for a — 0.5, i = 2 sec and uq = 0.1 cm. 



Mass (m) 
in a.m.u. 


Zpeak 

in cm 


(z) 
in cm 


30 


38.657865898827980 


38.65299987691765 


60 


38.657865903136155 


38.65299987691765 


90 


38.657865903934166 


38.65299987691765 


120 


38.657865904213070 


38.65299987691765 


150 


38.657865904342295 


38.65299987691765 



Here, i = p/m and p = —mg are obtained from Hamilton's equation of motion given by 

. _ dH 
dp ' 
dH 

dz 
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where, H = p^ /2m + mgz is the Hamiltonian of the system. 

The corresponding position distribution function is given by 

1*00 
Pc{z,t) = / D^{z,p,t)dp. 



So, the expression for the average value of z using Wigncr distribution function is given by 

, \ !T^zpc{z,t)dz 



lZoPciz,t)dz 
2 + a2 fl-e-T^ 



1 2 1 2 

ut - -gt ^ zo + ut- -gt 



which is exactly equal to the quantum mechanical expectation value of z. 
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FIG. 3: The variation of probability Pi{m) with mass (in a.m.u.) for different constant values of a. We take u = 10'^ cm s ^ , 
CTo ~ 10~ cm and e = ctq- 

The particles are projected upwards against gravity. Cassicaly, the refrence point z = reaches the maximum 
height z = Zjnax — uti — \gt\ at the time t = ti = u/ g. So, the probability of finding the particles Pi(m) within a 
very narrow detector region {— eto + e) around Zmax is given by 



Pi(m) 



r+c 
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(18) 



At a later time 1 = 1-2 = 2u/5, the classical peak of the wave packet returns to it's initial projection point at z = 0. 
Similarly, the probability of finding the particles P2{m) within a very narrow detector region {— eto + e) around 
z = is given by 



P2(m) 



\ij{z,t2)\^dz. 



(19) 



From Fig. |3]and Fig. 21 it is clear that both the probabilities arc mass dependent for an initial non-Gaussian position 
distribution for smaller masses. This effect of mass dependence of the probabilities occur essentially as the spreading 
of the wave packet under linear gravitational potential depends on the mass of the particles. These probabilities 
become saturated in the limit of large mass. As a increases, both the probabilities decrease for a particular value 
of 771. The departure from the result for the Gaussian (a = 0) case increases with the value of the non-Gaussian 
parameter, signifying a stronger violation of the weak equivalence principle. 



III. MEAN ARRIVAL TIME 



Now, Let us consider the quantum particles subjected to free fall from the initial position at z 
with the initial state given by Eq.(l) and with m = 0. Using the probability current approach [1^ 
time of the quantum particles to reach a detector location at z = Z is given by 



T = 



^^\J{z = Z,t)\tdt 

!^\j{z = z,t)\dt 
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4: The variation of probability P2{m) with mass (in a.m.u.) for different constant values of a. We take u - 



10 cm and e : 



: CTO- 



where, J{z,t) is the quantum probability current density. However, we emphasize that the definition of the mean 
arrival time used in the above equation is not a uniquely derivable result within standard quantum mechanics. It 
should also be noted that J{z,t) can be negative, so we take the modulus sign in order to use the above definition. 
From Eq. ((20)) , it can be seen that the integral of the denominator converges whereas the integral of the numerator 
diverges formally. To avoid this problem, several techniques have been employed [19| . In this paper, to obtain 
convergent results of the nu merator we use t he simple technique of getting a cut-off at i = T in the upper limit of 
the time integral with T = ^J2(\Z\ + iaT)/g- ot is the width of the wave packet at time T. Thus, our calculations 
of the arrival time are valid upto the Sctt level of spread in the wave function. 



Here, the expression for the Schrodinger probability current density J{z,t) with the initial non-Gaussian position 
distribution is given by 



J(z,t) 






'' ( 

— r V 



(0) 



arj 



(1) 



aV' 



(21) 



where 

77(0) = 8e 
•,(1) 



£2 



n^t 



-or 
2^ 



ry-' = r;(°)e"'^^ (e^^sinAi 



+ Am CTg [u — gt) , 

e^-" sinAz) - 2TTh^tao (e^'' cosAi 4- e^'' cos ^2) 



AirmTiag [e^"" sinAi - e^-" sinA2) 



V 



(2) _ 



2h^t ( z gf ) + Sto^CTq (u - gt) 



(1 



.S5 



Ee, 



- 2e^« COS A3) - 27r?i^tcroe " sin A3 -t- 2nmhc7'^{-l + e 



„B5^ 



and where the iS^'s are defined below Eq.(7). 

It is clear that J{z, t) is explicitly mass dependent. So, f is also mass dependent. From Fig. [SJ we can observe that 
for smaller masses f is mass dependent and for larger masses it becomes mass independent. For a particular value of 
mass, f increases with a. 

From Fig. [6l we can see that for large masses, f becomes almost constant with a. For a particular value of a, f 
decreases with increasing value of mass m. 

It is clear that the mass-dependence of the mean arrival time gradually vanishes as the mass is increased. Thus, 
compatibility with the equivalence principle emerges in the limit of large mass, or classical limit. Note that with the 
increase of the parameter a, one needs to go to higher values of mass in order for the mass-dependence to become 
insignificant. Thus, a non-Gaussian wave packet can exhibit non-classical features in a mass range where its Gaussian 
counterpart behaves classically. 
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FIG. 5: The variation of moan arrival time with mass (in a.m.u.) for different constant values of a. We take ao — 10 ^ cm, 
Z — —1cm and it = 0. 
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FIG. 6: The variation of mean arrival time with a for different constant values of m (in a.m.u.). We take ao = 10 ® cm, 
Z = —1cm and u = 0. 



IV. BOHM TRAJECTORY APPROACH 



Within the Bohmian interpretation of quantum mechanics [IJ] each individual particle is assumed to have a definite 
position, irrespective of any measurement. The equation of motion of the particle is, 



mx ^ -V(y + Q)|5=s, 
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(22) 



where Q = — ^ — -^ is called 'quantum potential energy' and it is explicitly mass-dependent and x(t) determines 
the path of the particle. In the gravitational potential V = mgz the quantum version of Newton's second law is given 

by 
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FIG. 7: A selection of Bohmian trajectories for u = (free fall) and ao = lO^^cm. In the left column tunable parameter a has 
the fixed value a — 0.5 and in the right one mass is constant and equal to tti = 100 amu. The solid black trajectory starts at 
{z}{0) , the dashed blue one at (z)(0) — 2ao and the dotted red one at (z)(0) + 2ao. 

To quote from Holland |5|, "It is immediately clear that the assumption of the equality of the inertial and passive 
gravitational masses does not imply that all bodies fall with an equal acceleration in a given gravitational field, due 
to the intervention of the mass-dependent quantum force term." Thus, WEQ is violated. The arrival-time problem is 
unambiguously solved in the Bohmian mechanics, where for an arbitrary scattering potential V{x), one finds [20[ that 
for those particles that actually reach x = X, the arrival-time distribution is given by the modulus of the probability 
current density, i.e., \J{X,T)\. 

The quantum potential Q is nonlocal, and thus it is expected that its effect on the violation of WEQ is manifested 
in a way similar to the violation observed as a consequence of the spread of the wave packet. We compute a set of 
Bohmian trajectories corresponding to free fall of our non-Gaussian wave packet with the tunable parameter a. Fig. 
[7] shows a selection of Bohmian trajectories exhibiting their spread depending upon the values of m and a. Note that 
for small mass, the trajectories with initial positions on the left and right of the center of the wave packet (mean 
(z)) spread out with evolution, indicating violation of WEQ. The magnitude of spread increases with a, signifying 
stronger violation of WEQ with increased departure from Gaussianity. The spread of the trajectories decreases as 
mass is increased, leading to the emergence of WEQ in the classical limit. 



CONCLUSIONS 



To summarize, in this work we have studied the violation of the gravitational weak equivalence principle in quantum 
mechanics using a non-Gaussian wave packet. The wave packet is constructed in such a way, that its departure from 
Gaussianity is represented by a continuous and tunable parameter a. The asymmetry of the wave packet entails a 
difference between its mean and peak, and causes the peak to evolve differently from a classical point particle, whereas 
the mean evolves like a classical point particle modulo a constant depending upon the parameter a. Such a result 
is consistent with the Ehrenfcst therorcm, as we have shown, and re-confirmcd using the Wigncr distribution. The 



11 

correspondence with the rcsuhs foUowing from a Gaussian wave packet is achieved in the limit of a ^- 0. 

The weak equivalence principle is shown to be violated first through the dependence on mass of the position 
detection probabilities of freely falling particles. The magnitude of violation increases with the increase of the 
non-Gaussian parameter a, signifying stronger viotation of WEQ with larger deviation from Gaussianity. We next 
compute the arrival time of freely falling wave packets through the probability current distribution corresponding 
to the non-Gaussian wave packet. The violation of WEQ is again exhibited with the dependence on mass of the 
mean arrival time. It is observed that in the limit of large mass the classical value for the mean arrival time is 
approached, thereby indicating the emergence of the WEQ in the classical limit. In this context it is worthwhile to 
mention that though our work follows as a natural consequence of quantum mechanics or quantum theory in the 
non-relativistic limit, it does not imply a violation of general covariance in the relativistic domain. An interesting 
connection between the non-relativistic limit of quantum theory and the principle of equivalence has recently been 
discussed [21|. Finally, a computation of Bohmian trajectories further establishes the stronger violation of WEQ 
by a non-Gaussian wave packet. We conclude by noting that realization of exactly Gaussian wave packets is rather 
difficult in real experiments. One of the consequences of our results is to enable relaxation of the wave packet to be 
Gaussian, and therefore it can facilitate the experimental observation of the violation of WEQ, as well as enable a 
quantitative verification of the way the violation of WEQ depends on the departure from the Gaussian nature of the 
wave packet. Other ramifications of this result call for further study. 
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